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Abstract.  A  convex  body,  with  high  thermal  conductivity,  is  immersed  in  a  nonuniformly  flowing  gas.  The 
body  is  small  compared  to  the  mean  free  path,  which  in  turn  is  small  compared  to  the  macroscopic  length 
scale  of  the  gas.  The  equilibrium  temperature  Tw  of  the  body  is  calculated.  For  an  axially  symmetric  body 
in  a  simply  shearing  gas  of  temperature  T  one  obtains  the  equilibrium  temperature 

Tw  .  /fo  Pxy  .  2  /)  •  /n  \ 

—  =  1  +  —  — ^  sm  0sin(2¥)). 

1  op 

(this  is  for  the  case  that  the  body  is  at  rest  with  respect  to  the  gas).  0,  <p  are  polar  angles  of  the  axis  of  the 
body  ( z  is  the  polar  axis),  a  is  a  geometric  shape  factor  of  the  body  (which  vanishes  for  a  sphere)  and  fi 
depends  on  the  Sonine  coefficients.  /3  takes  the  value  1  if  only  the  lowest  order  Sonine  term  is  retained,  p 
is  the  pressure  and  pxy  (the  non-vanishing  component  of)  the  viscous  pressure  tensor. 


INTRODUCTION 

The  phenomenon  of  thermophoresis  -  the  force  acting  on  a  small  body  in  a  gas  with  a  temperature  gradient 
-  is  in  the  limit  of  freee  molecular  flow  caused  by  the  non-Maxwellian  nature  of  the  molecular  distribution 
function,  see  [1].  There  is  an  analogous  phenomenon:  in  a  shearing  gas,  due  to  the  non-Maxwellian  character 
of  the  molecular  distribution  function,  there  will  be  a  heat  transfer  to  a  small  body.  -  For  a  general  review  of 
free  molecular  flow,  see  Schaaf  [2].  See  also  Sone  [3]. 

Bell  and  Schaaf  [4] ,  see  also  Schaaf  [2]  calculated  the  heat  transfer  to  an  infinite  circular  cylinder  in  the  free 
molecular  flow  approximation.  They  used  the  first  order  Chapman-Enskog  distribution  function.  As  a  result 
they  found  a  contribution  to  the  heat  transfer  from  the  shearing  of  the  gas.  Consequently,  the  equilibrium 
temperature  of  the  body  was  also  shown  to  be  affected  by  the  shearing  of  the  gas. 

In  the  present  work,  we  consider  a  convex  body,  which  is  small  compared  to  the  mean  free  path.  The 
contribution  from  the  mean  flow  of  the  gas  is  calculated  up  to  second  order  in  the  Mach  number  and  the 
contribution  from  the  shearing  is  calculated  to  first  order  in  the  Mach  number  and  first  order  in  the  Knudsen 
number.  We  assume  free  molecular  flow,  see  Kogan  [5].  We  also  use  the  common,  test  body,  approximation  of 
neglecting  the  influence  of  the  body  on  the  incoming  flow  of  gas  molecules. 

We  use  the  Maxwell  model  for  the  surface  assuming  that  a  fraction  aT  of  the  incoming  particles  is  diffusively 
reflected,  whereas  the  fraction  (1  —  aT)  is  specularly  reflected.  We  further  assume  that  the  surface  of  the  body 
is  characterized  by  the  energy  accomodation  coefficient  ae.  See  Kogan  [5].  The  gas  is  taken  to  be  monatomic. 

The  equilibrium  temperature  is  first  derived  for  a  general  distribution  function.  This  is  first  applied  to  a 
homogeneous  flow  and  then  to  a  shearing  flow  with  zero  velocity  at  the  body.  The  case  of  a  shearing  flow  with 
non-zero  velocity  at  the  body  is  in  the  approximation  used  given  by  superposition  of  the  two  contributions. 


GENERAL  DISTRIBUTION  FUNCTION 
Net  influx  of  kinetic  energy 

Let  us  now  start  by  writing  down  some  relations  of  general  validity.  See  Kogan  [5] . 
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The  Maxwellian  distribution  is  written 


Here,  n  is  the  number  density  of  molecules,  kB  the  Boltzmann  constant,  T  the  temperature,  and 
molecular  mass.  Ci  is  the  dimensionless  velocity 


Ci  = 


m 

2  kBT 


Ci, 


where  C,  is  the  ordinary  molecular  velocity. 

For  a  general  distribution,  the  total  influx  of  particles  is 


N(i)  =  -[  (C  ■  n)/d3C. 
J  n-C<0 


In  particular,  for  a  Maxwellian  distribution  N ^  is  given  by 


Nm  =  n 


kBT 

2ivm' 


The  conservation  of  particles  at  the  surface  is  expressed  by 


N(l)  =  nr  \  I  ^bTt  =  n„ 


2tt  m 


kBTw 

2irm 


Index  r  refers  to  reflected  and  w  to  wall. 
The  total  influx  of  kinetic  energy  is 


EW  =  ~  /  C2(C.n)/d3C. 

^  J n-C<0 


In  particular,  for  a  Maxwellian  it  takes  the  value 

E  =  rJ^_  (2 kBT  3/2 
2  '  irm 

According  to  the  Maxwell  model,  the  outflux  of  particles  is  given  by 

E(r)  =  (1  -  aT)E(i )  +  aTEM{Tr). 
Thus  the  net  influx  of  kinetic  energy  is 

E  =  E(i)  —  i?(r)  =  otr(E(i )  -  Em{Tt)). 
Let  us  also  introduce  the  correspondence  of  (8)  for  wall  conditions 

E(w)  ~  (1  ®'r)E(i)  +  OItEm{Tw). 

The  energy  accomodation  coefficient  is  given  by 

E(i)  E(r) 


E(i)  E(.w) 


Using  (5)  we  find  the  net  influx  of  kinetic  energy  as 


aeaT{E(i) 


E  =  aeaT(E(i)  —  Em{Tw)) 


Tw  nrTr/2 
T  nT1/2 


Em)  =  OLeaTEM{ 


Ed) 

Em 


Tw  IV(i) . 
T  Nm> 


(1) 

m  the 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 
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Equilibrium  temperature  for  high  conductivity  body 


Let  us  now  for  simplicity  assume  that  ae  and  aT  are  constant.  We  further  assume  that  the  thermal  conduc¬ 
tivity  of  the  body  is  so  high,  that  Tw  is  uniform  in  the  body.  The  total  influx  of  kinetic  energy  to  the  body 
can  then  be  written  ( S  is  the  area  of  the  body) 

J Eds  =  aeaTEMS[( |g)  -  (13) 

The  overbar  denotes  surface  average 

7  =  |//d, 

At  stationary  conditions,  the  (equilibrium)  temperature  of  the  body  is  thus 


Tw  _  . 

T  KEM)nNM> 


(14) 


HOMOGENEOUS  FLOW 


Let  us  now  consider  a  homogeneous  flow  of  the  gas.  See  Bird  [8].  We  introduce  the  dimensionless  mean 
velocity  (e  is  a  unit  vector) 


Vi  =  Vei  = 


m 

2 k^T 


Vi. 


The  speed  of  sound  is  given  by 


C.s  — 


'jksT 


m 


Here,  7  =  5/3  for  a  monatomic  gas.  The  Mach  number  is  then 

Ma=y|v. 

The  distibution  function  is  a  Maxwellian,  where  Ci  is  replaced  by  Ci  —  V,;. 

The  particle  influx  is  now  calculated  as 

N(i)  =  AriW[exp(— V/)  -  v/7rV„(l  +erf(-V„))]. 

If  we  restrict  ourselves  to  terms  up  to  second  order  in  the  Mach  number,  we  can  use 


(15) 


(16) 


(17) 


(18) 


erf(a;)  «  — =x, 

-nA 


giving 


^0 

Nm 


=  l  +  V2n-^Vn. 


We  find  for  the  surface  average  (the  term  linear  in  V  gives  no  contribution  as  the  surface  is  closed) 


(^)  =  l+V2ATijeiej. 

Here,  the  purely  geometric  tensor  Nij  is  ( S  is  the  area  of  the  body) 


(19) 
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Nu  =  —  riiTiids. 


It  is  convenient  to  split  Nij  into  its  irreducible  components  (<  ...  >  denotes  symmetric  traceless  part) 


The  eigenvalues  of  lie  in  the  interval  [—1/3, 2/3].  Let  us  note  that  N<ij>  vanishes  for  a  sphere. 

For  an  axially  symmetric  body,  with  axis  N,  the  tensor  N<ij>  is  given  by  a  scalar  a, 

N<ij>  = 


3  1 
a=2  S 


/[( n-N)2-i]ds 


For  a  flat  body  a  =  1.  For  a  needle  shaped  body  a  =  — 1/2. 
Similarly  we  have 


E{i)  =  Em{{~z  hi)  exp(— V^)  -  \ftVn(^r  +  |)[1  +  erf(-V„)]}. 


Approximately, 


The  surface  average  is  thus 


2  4' 


E(i)  V2  3  ,  5 

E^=1  +  Y  +  2Vn~4^Vn- 


(pr^)  —  1  +  V2(l+-AT<jj>eiej). 


Inserting  (19)  and  (25)  into  (13)  we  find  the  total  influx  of  kinetic  energy  for  a  high  conductivity  body.  (14) 
gives  the  equilibrium  temperature  to  second  order  in  the  Mach  number  as 

Tw  5Ma2  .  3 .  .  . 

■jT  =  H - g (1  +  ^<ij>eiej)-  (26) 


Axially  symmetric  body 

Let  us  consider  an  axially  symmetric  body.  If  there  is  no  external  heat  transfer  to  the  body,  in  stationary 
conditions,  it  will  have  the  temperature  (with  an  error  V3)  -  d  is  the  angle  between  the  axis  of  the  body  and 
the  direction  of  the  flow  - 

~fr  =  1  +  — g — [1  +  ^(Scos2  0  -  1)].  (27) 


For  a  sphere, 


^  =  1  +  iMa2. 

T  9 


This  is  also  the  mean  temperature  if  a  small  body  of  arbitrary  shape  is  randomly  oriented. 

For  a  needle  shaped  body 

=  1  +  \Ma2(^  -  \  cos2  0)  (29) 
In  particular,  if  the  axis  of  the  body  is  orthogonal  to  the  flow, 

^  =  1  +  \Ma2.  (30) 

For  an  infinite  cylinder  with  axis  orthogonal  to  the  flow,  this  coincides  with  the  result  of  Bell  and  Schaaf  [4]. 
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SHEARING  FLOW 


Chapman-Enskog  solution 

In  a  pure  shearing  flow,  according  to  the  Chapman-Enskog  method,  the  distribution  function  is  given  by,  see 
Chapman  &  Cowling  [6] 


/  =  /«»(  i  +  4>). 


4>  —  ——B(C2)C<iCj>v<ij>  . 


(31) 


B  is  usually  expanded  in  Sonine  polynomials 


B(C2)  =  Y/bnSif2(C2). 


(32) 


n=0 


The  coefficient  6q  is  related  to  the  viscosity, 


/x  =  ksTbo 


Let  us  introduce 

B(C2)  =  B(C2)/b0 


Using  the  gas  law  p  =  nksT,  we  can  write  ( p<ij>  =  —2p,v<ij>  is  the  viscous  pressure  tensor) 

<p  =  -B(C2)C<iCj>p<ij>  . 


(33) 


Particle  flux 

The  shearing  gives  the  following  contribution  to  the  influx  of  particles 

p  J 

Here,  we  have  introduced  the  dimensionless  integral 

Iij=—  [  CnB(C2)C<iCj>e~c2  d3C  =  lNn<irij> 
n  Jc„<o 

The  last  expression  follows  from  symmetry.  Contracting  twice  with  the  normal  vector  and  introducing  spherical 
polar  coordinates,  we  find,  writing 

Cn  =  Ck-rik,  (34) 

In  =  — —  f  Cn{C2n  -  l-C2)B{C2)e~c2 d3C 

7r  Jcn<o 

=  --pNK.  (35) 

7T 

Here, 

/»OC 

Pn  =  /  B(C2)C5e~c2dC  (36) 

Jo 

1  - 

=  -  I  B(x)x2e~xdx. 

2  Jo 
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Note  that  /3jv  =  1,  if  only  the  first  term  of  the  expansion  of  4>  in  Sonine  polynomials  is  retained.  This  is  usually 
a  good  approximation.  -  Further,  the  angular  integral  is 

/*7T  1 

K  =  2tt  I  (cos2  0  —  —)  sin  0  cos  0d6  =  —  —  . 

J  tt/2  3  6 


We  conclude  that 


In  =  2&n- 

The  contribution  from  the  shearing  to  the  inflow  rate  of  particles  is  thus 

SKMi>n<ini>NM. 

A  p 

The  total  inflow  of  particles  is  then  given  by 

N(i)  =  Nm(  1  +  ^E^n<inj>) 

z  p 


Hence, 


(37) 


(38) 


Flux  of  kinetic  energy 

Let  us  write  the  contribution  from  the  shearing  to 


1  P<ij>  TEi 


2  p 


Here,  the  tensor  7®  is  (the  expression  to  the  right  follows  from  symmetry) 

iE  =  ~  [  C2CnB(C2)C<iCj>e~c2 d.3C 

n  Jcn< o 


Contracting  Iij  with  n^rij  we  find,  using  spherical  polar  coordinates, 

Ie  =  -JeK. 

7 T 


In  this  expression,  the  radial  velocity  integral  is 


Je 


Let  us  introduce  a  dimensionless  coefficient 

1  r°°  ^ 

/?£  =  -/  B(x)x3e~xdx.  (39) 

6  Jo 

Usually,  the  n  =  0  term  in  (32)  is  the  dominating  one.  If  the  other  terms  are  neglected,  @e  takes  the  value  1  . 
-  Je  can  now  be  written 


Je  =  3  (3e- 
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Hence 


Thus  the  total  rate  at  which  heat  is  incoming  is 

Ed)  =  (!+  ^ EP<l^>n<inj>)EM  (40) 

We  can  easily  estimate  the  size  of  the  contribution  from  shearing.  We  have  /i/p  ~  X/cs,  where  A  is  the  mean 
free  path  and  cs  the  speed  of  sound.  If  we  estimate  u<jj>  ~  v/l,  we  find  the  relative  size  of  the  contribution 
to  be  Ma  Kn,  ( Ma  =  v/cs  the  Mach  number  and  Kn  =  A /l  the  Knudsen  number.) 

The  total  incident  flux  is  thus 


(|g>  =  (l  +  ^^V<«>).  (41) 

Total  rate  of  heat  transferred  to  the  body 

Body  of  arbitrary  shape 

The  total  influx  of  kinetic  energy  is  now  given  by  (13). 

We  see  that  also  for  Tw  =  T,  there  is  a  net  influx  given  by 

^aeaT-EMrc<»rv>P<p  >  ■  (42) 

We  have  here  introduced 


(3  =  3  (3b  —  2/3  a'  ■ 


(43) 


Note  that  /3  =  1,  if  only  the  first  term  in  the  Sonine  polymial  expansion  of  0  is  retained. 

Let  us  now  consider  a  body,  which  adjusts  its  temperature  so  that  the  total  inflow  vanishes.  To  a  first 
approximation  in  the  shearing, 


T 

±W 

~T 


=  1  +  yP<o> 

4  p 


N<ij>7 


(44) 


where  /3  is  given  by  (43). 


Axially  symmetric  bodies 

Let  us  in  particular  consider  an  axially  symmetric  body  in  plane  shear  flow.  In  suitable  coordinates,  the 
only  non- vanishing  component  of  the  velocity  gradient  is  vx>y .  Let  us  denote  the  polar  angles  (z  is  the  polar 
axis)  of  the  axial  direction  of  the  body  by  6 ,  ip.  Then 

N<ij>P<ij>  =  ^PxV  sin2  6  sin(2<p). 

In  particular,  for  a  body  free  to  adjust  its  temperature  until  the  total  heat  flow  vanishes,  the  temperature 
of  the  body  relaxes  to 


=  1  +  ^  —  sin2  0sin(2p).  (45) 

i  op 

For  a  needle-shaped  body,  putting  /3  =  1  (neglecting  higher  Sonine  polynomials),  this  coincides  with  the  result 
of  Bell  and  Schaaf  [4], 
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A  needle-shaped  body  will  then  settle  down  to  a  temperature,  which  is  higher  than  that  of  the  surrounding 
gas,  if  the  projection  onto  the  plane  of  shearing  of  the  axis  of  the  body  is  in  the  first  quadrant,  but  at  a 
temperature,  which  is  lower  than  that  of  the  gas,  if  the  axis  is  in  the  second  quadrant. 

Let  us,  finally,  consider  a  body,  which  is  kept  at  the  temperature  of  the  surrounding  gas.  We  then  have  the 
total  rate  at  which  heat  is  transferred  to  the  body 

^~aeaTEMS^^-  sin2  $sin(2<^).  (46) 

We  conclude  that  needle-shaped  bodies  absorb  heat,  when  the  projection  of  their  axes  on  the  plane  of  shearing 
lie  in  the  first  quadrant.  They  emit  heat  when  their  axes  are  in  the  second  quadrant.  For  flat  bodies  the 
situation  is  reversed. 


TRANSLATION  AND  SHEARING 


For  an  arbitrary  flow  field  with  shear,  we  can  to  the  order  considered  add  the  two  contributions  to  —  1 
to  find  (we  recall  that  e  is  the  unit  vector  in  the  direction  of  the  flow) 


-TjT  —  1  +  g  Afa2  +  (— Ma?eiej  +  )N<ij>. 


(47) 


This  is  with  an  error  Ma3,  KnMa2,Kn 2  Ma. 

Let  now  also  consider  an  axially  symmetric  body  in  a  flow  field,  where  the  only  non-vanishing  velocity 
component  vx  is  a  function  of  y.  This  means  that  e  •  N  =  cosd  is  replaced  by  sin  0  cos  ip.  We  obtain 


Tw  5Ma2  r  a  f  .  ^  n  2  m 

—  =  !  + -y-[!+ -(3sin  6>cos  ip-1)] 


Totally, 


r^'w  1  i  'iMci  .  d .  .  2  nr 

-T  =1  +  — (1-4)  +  a,m  81  12 


2«'—  cosV+F™™(V)]. 


Ppx 

8  p 


(48) 
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